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Rachunek zaburzeń Hartree-Focka w skończenie wymiarowych bazach funkcyjnych. 
Część I. Formalizm rachunku zaburzeń w ustałonych i zmiennych bazach funkcyjnych 
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The Hartree-Fock Perturbation Theory for Finite Basis Sets. 
Part I. Formalism of Perturbation Theory for Fixed and Variable Basis Functions 


WSTĘP 


Ogólnie wiadomo, że tylko nieliczne problemy kwantowo-mechaniczne mogą być roz- 
wiązane w sposób ścisły [1, 2]. Już obliczenie funkcji falowych dla izolowanych atomów 
wieloelektronowych i molekuł napotyka na szereg podstawowych trudności i wymaga 
stosowania metod przybliżonych [3]. Jeszcze bardziej kłopotliwe jest uzyskanie informa- 
cji o wpływie na te układy czynników zewnętrznych, takich jak pola elektryczne i mag- 
netyczne [4]. Warto podkreślić, że większość naszych wiadomości o atomach i moleku- 
łach pochodzi właśnie z badań ich oddziaływania z polami zewnętrznymi. 

Źródłem informacji o wpływie zewnętrznego pola, reprezentowanego przez operator 
zaburzenia AHO), na dany układ wieloelektronowy opisywany hamiltonianem H), są 
rozwiązania równania 


[A +AH WA) = EA) VA) (1) 
Parametr A jest traktowany jako miara wielkosci fizycznego pola zewnetrznego (lub 
.sktadowej pola) i jest jednoczesnie parametrem rozwiniecia w postepowaniu perturba- 


cyjnym [5]. Traktujac zaburzenie MU) jako mate [1, 2] oraz zakładając znajomość 
dokładnych rozwiązań W(9 i E(0) równania Schrödingera dla układu niezaburzonego 


HO) yo) = pio) (0) 
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można rozwiązania równania (1) aproksymować przez szeregi potęgowe względem A [1,6] 


W = YO) LAVO + 2G) +... 
E = BO EO REO +... (2) 


Liczne własności układów wieloelektronowych, charakteryzujące ich oddziaływanie 
z zewnętrznymi połami, proporcjonalne są do członów rozwinięcia enezgii, czyli do po - 
prawek k-tego rzędu EG), Zależnie od charakteru zaburzenia (np. pole elektryczne), mogą 
to być multipolowe (dipolowe, kwadr-upolowe, ...) momenty p, polaryzowalnošci a 
i hiperpolaryzowalnošci 6 [7] 


KE), a- EO), pe EO), 

Wymienione tutaj wielkości w pełni tłumaczą duże zainteresowanie metodami pertur- 
bacyjnymi, które dają możliwość ich dokładnego obliczenia. 

Do obliczenia poprawki pierwszego rzędu Æ (1) wymagana jest jedynie znajomość 
funkcji falowej rozpatrywanego stanu układu niezaburzonego [1, 6]. Problem obliczenia 
wyższych poprawek energetycznych związany jest z uprzednim wyznaczeniem odpowied- 
nich poprawek do funkcji falowej [1, 6]. Te ostatnie są rozwiązaniami niejednorodnych 
równań różniczkowych, które otrzymuje się po podstawieniu równań (2) do równania (1) 
i separacji według potęg parametru A. Bezpośrednie, dokładne rozwiązanie tych równań 
możliwe jest tylko w nielicznych przypadkach o niewielkim znaczeniudlatypowych probie: 
mów chemicznych [8]. W ogólnym przypadku można poszukiwane poprawki Y( ) roz- 
winąć w szereg funkcji własnych operatora H(9), Prowadzi to do tzw. wyrażeń spektral- 
nych na poprawki do funkcji falowej i do energii [8]. W tych formalnie prostych wyraże- 
niach występuje problem sumowania po pełnym widmie stanów własnych operatora 
HO), Ponieważ w praktyce nie dysponujemy pełnym widmem operatora HO), stąd też 
zastosowanie rachunku zaburzeń w tym standardowym ujęciu jest bardzo kłopotliwe. 

Bezpośrednie sumowanie po układzie zupełnym stanów własnych może być wyelimi- 
nowane przez wprowadzenie sformułowań wariacyjnych dla odpowiednich równań ra- 
chunku zaburzeń [9]. Przejście do sformułowania wariacyjno-perturbacyjnego nie ozna- 
cza jednak usunięcia trudności numerycznych. W praktyce, chcąc uzyskać na podstawie 
teorii informacje liczbowe o własnościach układów wieloelektronowych, naležy wprowa- 
dzić pewne założenia upraszczające. Jednym z nich jest zastąpienie dokładnych rozwią- 
zań Y) j EE) przez rozwiązania wynikające z zastosowania przybliżenia Hartree-Focka. 


RACHUNEK ZABURZEŃ W PRZYBLIŻENIU HARTREE-FOCKA 


Omówienie rachunku zaburzeń Hartree-Focka ograniczone zostanie do pierwszego 
rzędu względem funkcji falowej. Podane tutaj sformułowanie zaczerpnięto z prac [12, 13]. 

Rozpatrzmy 2n — elektronowy układ zamknięto-powłokowy w stanie podstawowym, 
ppisywany w przybliżeniu Hartree-Focka pojedynczym wyznacznikiem Slatera skonstru- 
awanym z n orbitali uf) spełniających równania Hartree-Focka 


SA A A 7 
KO) ON) = HOMO 
Hamiltonian h©)(1) dla układu niezaburzonego określony jest wyrażeniem 
HO) = HO) + 5 tu; IG 77 lu (2) 


w którym f(2)(1) jest suma operatorów energii kinetycznej i energii potencjalnej elektro- 
nu w polu jąder 


fONM=>WVI + MY) 


(Gia = R2(2—P,). ; 
Operator Pj, jest operatorem permutacji elektronów 1 i 2, a indeks 2 w wyrażeniu (+++), 
oznacza całkowanie po współrzędnych elektronu 2. 
Jeżeli układ poddany zostanie zaburzeniu będącemu sumą zaburzeń jednoelektro- 
nowych À : 
i pat i 
AH) => OW, (3) 
t i=1 
to zaburzone zagadnienie własne w przybliżeniu Hartree-Focka dane będzie równaniem 
ACL; Nul; A) = er Ou; (1;X), (4) 
w ktörym 
WC) FD AND +Ż (IC i ON 
=1 
Rozwijając u; (X), e;(X) i RAJ w szeregi potegowe względem A i podstawiając je do rów- 
nania (4), a następnie dokonując separacji rzędów względem potęg parametru A, otrzy- 
mamy równania Hartee-Focka w kolejnych rzędach rachunku zaburzeń. W pierwszym 
rzędzie odpowiednie równanie ma postać 
[AO — e; Ou; = [e;0) — 01) mg DO), (5) 
gdzie 
e) = (4; (O |) jy) 
n 
gO =,2, [lu [Gu luf?» + Gul 1617 14), ]. 


Rozwiązanie równania (5) daje funkcję u;'') zwane najczęściej poprawkami pierwszego 
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rzędu do orbitali Hartree-Focka. Znajomość tych funkcji wystarcza do obliczenia wiel- 
kości E(?), która oznacza teraz poprawkę drugiego rzędu do energii całkowitej w przyb- 
liżeniu Hartree-Focka [14] 


E = 2 590, „to XE (uz OYO uz), 

Rozwiazanie röwnania pierwszego 1zedu O wymaga znajomosci poszukiwanych 
funkcji O), co wynika z postaci operatora g! . Tak więc mamy do czynienia w rzeczy- 
wistości ze sprzężonym układem równań i stąd omawiana metoda nosi nazwę sprzężone- 
go schematu rachuku zaburzeń Hartree-Focka (CHF = Coupled Hartree-Fock perturba- 
tion theory) [14, 15]. 

` Najprostszą metodą formalnego rozwiązania równania (5) jest rozwinięcie eee 
nych funkcji u; jU) na układ zupełny UT ha orbitali Hartree-Focka u; (0) Zato- 
zenie znajomošci ukladu zupelnego orbitali u; (0) jest jednak praktycznie niemožliwe do 
speinienia. Przyjmujac, Ze röwnania Hartree-Focka dla problemu niezaburzonego sa 
rozwiązane ściśle dla n zajętych orbitali uj(?) układu, można poszukiwać rozwiązań rów- 
nania (5) przez wprowadzenie zasady wariacyjnej. Prowadzi to do wyznaczania funkcji 
uj) poprzez minimalizacje odpowiedniego funkcjonatu, dla ktörego röwnanie (5) jest 
röwnaniem Eulera [16] 


TO) (450) = (uy (OI HOO — e O bug) + Xu OYO — e; O bj) + 
+ (2-87) Da O 16 17 kag Pa O + (up Oj IG lu; uj OYJ. 
=1 2 


W konkretnych obliczeniach metoda CHF, zaröwno w wariancie perturbacyjnym, jak 
i wariacyjno-perturbacyjnym, zmuszeni jestesmy do zastapienia doktadnych orbitali 
Hariree-Focka orbitalami przyblizonymi, uzyskanymi metoda pola samouzgodnionego 
(SCF= Self-Consistent Field) w skoriczenie wymiarowej bazie znanych funkcji. W pertur- 
bacyjnym schemacie CHF (5) poszukiwane funkcje u;“) rozwijane są na skończony 
układ przybliżonych orbitali SCF, Schemat wariacy INY stwarza natomiast 
większe możliwości wyboru w postaci funkcji u;‘*’ [12], co jest istotną jego zaletą. 

Zastąpienie schematu perturbacyjnego przez wariacyjno-perturbacyjny nie usuwa 
jednak podstawowych trudności metody CHF. Wynikają one z konieczności stosowania 
techniki iteracyjnej oraz obliczania znacznej liczby całek dwuelekronowych [12]. Ograni- 
cza to w pewnym stopniu stosowalność metody w obliczeniach molekularnych do stosun- 
kowo niewielkich układów [4, 17]. W obliczeniach dla układów większych wprowadza 
się niekiedy uproszczenia, które polegają na częściowym lub całkowitym zaniedbaniu 
członów sprzęgających w funkcjonale CHF. Powoduje to jednocześnie częściową lub 
całkowitą eliminację całek dwuelektronowych. Zależnie od stopnia zaniedbania członów 
sprzęgających, w równaniach metody CHF otrzymuje się różne metody zwane rozsprzężo- 
nymi schematami rachunku zaburzeń Hartree-Focka (UCHF = Uncoupled Hartree-Fock 
perturbation theory)[16]. 

Przejście do przybliżonych rozwiązań SCF, wykorzystujących skończenie wymiarowe 
bazy funkcyjne, powoduje utratę dokładności wyników metody CHF [12]. Uzyskanie 
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wartości E (2) bliskich granicy Hartree-Focka dla energii drugiego rzędu wymaga stosowa- 
nia baz funkcyjnych o bardzo dużym wymiarze [18], co stwarza określone trudności na- 
tury numerycznej. Warto w tym miejscu dodać, że dokładność wyników E (2) uzyskiwa- 
nych metodami typu UCHF jest z reguły niższa od dokładności wyników CHF [19], co 
jest oczywiście konsekwencją dokonanych u park Tak więc, można ogólnie stwier- 
dzić, że uzyskanie dokładnych warości E 2) jest praktycznie niemożliwe metodami 
UCHF oraz bardzo kłopotliwe metodą CHF. 

Konieczność stosowania w obliczeniach CHF baz funkcyjnych o dużych wymiarach 
jest wynikiem powolnej zbieżności uzyskiwanych wartości E (2) do dokładnej energii 
drugiego rzędu Hartree-Focka [12, 14]. Ponadta wraz ze wzrostem wymiaru bazy rosną 
trudności natury numerycznej, wynikające między innymi z konieczności obliczania bar- 
dzo dużej liczby całek dwuelektronowych. 

Zbieżność w metodzie CHF (SCF) może być istotnie poprawiona, jeżeli w obliczeniach 
wykorzystane zostaną bazy nie ustalonych, ale zmiennych funkcji, to jest takich, które 
bezpośrednio zależą od wielkości zewnętrznego zaburzenia [x; = x;(A) J. 

Dysponując m wymiarową zmienną bazą funkcyjną, można przedstawić zależne od 
zaburzenia orbitale Hartree-Focka u;(1;A) = u;(X), wynikające z równania (4) w postaci 


4:09 = E ey Qs 
która pozwala zapisać u;() i uj jako 


uj = Z Oy, (0) 
ja 


uj = 3 er Oy (OO + p, (0) (1) 
2 lo Xj Be aaa 


gdzie x, 1) = =[(0x,A)/@A)]}, 0° 


Jak widać, przy zastosowaniu m wymiarowej, zmiennej bazy funkcyjnej, orbital ` 
Hartree-Focka pierwszego rzędu u; 1) rozwijany jest w Rate sa na układ 2m 
funkcji. Przechodząc natomiast do przypadku ustalonej bazy mamy X; ©) = Qi orbital 
pierwszego rzedu jest kombinacja liniowa tylko m funkcji X; (0) .Wobec | powyzszego moz- 
liwość poprawy zbieżności metody CHF poprzez wykorzystanie zmiennych baz funkcyj- 
nych jest rzeczywiście realna, pod warunkiem jednak, że funkcje [x;(A)] przy braku zew- 
nętrznego pola AS 0) beda stanowity dobra (z punktu widzenia doktadnej energii 
Hartree-Focka E(?)) bazę dla problemu niezaburz-onego. Ponadto efektywność zastoso- 
wania zmiennych baz w obliczeniach CHF zależy od charakteru zależności funkcji x;(A) 
od wielkości zaburzenia A TZ Będzie to dokładnie omówione w drugiej części tej 
pracy. 

W przypadku zmiennych baz funkcyjnych zachodzi również potrzeba uogólnienia 
schematu metody CHF [11]. 
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METODA CHF W ZMIENNYCH BAZACH FUNKCYJNYCH 


Metoda CHF (SCF) w ustalonych bazach funkcyjnych może być sformułowana bądź 
w tzw. ujęciu orbitalnym, które było przedstawione poprzednio, bądź też w ujęciu 
macierzowym [24]. Podejście orbitaine ma mariin sakay interpretacyjne, ale wyma- 
ga bezpośrednich obliczeń perturbowanych orbitali uj“). Ten etap nie wystę puje w uję- 
ciu macierzowym, co jest korzystne z numerycznego punktu widzenia. 

W tym wygodnym ujęciu macierzowym metoda CHF (SCF) została niedawno opraco- 
wana dla przypadku nieortogonalnych i bezpośrednio zależnych od zaburzenia baz funkcyj- 


nych [11]. 

Rozważany jest, jak poprzednio, układ zamkniętopowłokowy z n podwójnie obsa- 
dzonymi orbitalami [41, ... , 4n] = u. Orbitale te przedstawione są jako liniowe kombi- 
nacje m zmiennych, rzeczywistych funkcji bazy [X; , ..., Xm] = XX = xO) 

> u=xT | 
gdzie kolumny macierzy T(mxn) są wektorami własnymi macierzy operatora Hartee- 


-Focka i wynikają z rozwiązania macierzowego równania SCF HF przy spełnieniu warun- 
ków ortonormalności 


hT = STe 
T'ST =1. 


Diagonalna macierz € zawiera energie orbitalne, a S jest macierzą całek nakładania orbita- 
li bazy Siz = (x; by). Macierz T* jest transponowaną i sprzężoną macierzą T. 
Wprowadzajac macierz gestošci Focka-Diraca 
R=TT* 
można równania SCF i warunki ortonormalności zapisać równoważnie jako 


hRS = SRh, 
RSR=R. (6) 


Reprezentacja macierzowa operatora Hartree-Focka dana jest przez 
h =f +G(R.g), 
gdzie f reprezentuje w bazie y jednoclektronowy operator anap 
| f= f OO + gO), 


Macierz G reprezentuje potencjat Hartree-Focka i zalezy od macierzy gestošci R i men 
macieszy g zawierającej całki dwuelektronowe* 
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GR g); = 3 Revi ki 
s ij, kl = 26GX4 lxi Xp — KKK da XP. 
Catkowita energia elektronowa uktadu okreslona jest wyrazeniem 
E=2Tr[f+ 0,5G]R. | (7) 


Pamietajac, że każda funkcja bazy zależy explicite od wielkości zewnętrznego zaburzenia, 
można elementy wszystkich macierzy rozwinąć w szeregi potęgowe względem A 


S=S©) 4180) 47280) +... 
f= £0) 4) + 72000 4. © . . 
h=h® + Ah) + XA7h02) + ,,, (8) 
R5RO ARO + XRO +... 
G=GRO g) + AGRO, g0) + GRO g+... 


gdzie dla przykładu mamy w pierwszym rzędzie 


Sy M = 04V fy) + gal) by (1) 
fy? = = die) + Paos If! xj 2 AKO [so hy) 


ianalogicznie dla pozostatych wielkošci. W powyžszych wyrażeniach | 


x = x;(0) 
xf = [(0x; AVON), 0 


oraz görny indeks zero oznacza wielkosci dla problemu niezaburzonego. 

Podstawiając równania (8) do równań (6) otrzymamy po separacji rzędów względem - 
potęg A równania SCF z warunkami ortonormalności dla kolejnych vy rachunku za- 
burzen. Dla pierwszego rzędu równania te maja postać 


KOR@S@) + pORA)SO) + pOROMG® — SOROKO + SORO + SOROR 
ROS(ORO) + ROSOR) + ROSOR) = RO) 


Znając rozwiązanie dla problemu niezaburzonego RO), można z ostatnich równań wyz- 
naczy€ macierz gęstości w pierwszym rzędzie RO, Naležy jednak zauważyć, że macierz 
PR Hartree-Focka w pierwszym rzędzie zależy poprzez G od poszukiwanej macie- 
rzy RU) . Pociąga to za sobą konieczność iteracyjnego rozwiązywania równań lakid 
rzedu [11]. 

Podstawiajac rõwnania (8) do wyraženia D, otrzymamy dość złożone wyrażenia, 
określające kolejne poprawki energetyczne El . Wykorzystując wygodną i elegancką 
technikę operatorów rzutowych [11, 24] można wyżóńie te uprościć, otrzymując osta- 
tecznie dla E) j EO) 
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ED =2Ty[f") + 0,5G(R) g) -h ROSON RO) 


E) = 2 Ty[f£(2) n AO ROS) EE 0,5G(R(?) g? » me 0,5h! pg (0)g(1 Ro e 
+ Tr{ ft —pORMg@) + GRO g) —g(t RAC] RG A 


W analogiczny sposób można wyprowadzić równania określające wyższe poprawki 
energetyczne, przy czym z powodu załeżności funkcji bazy od wielkości zaburzenia będą 
one coraz bardziej złożone. Na pie w wyrażeniu na Æ (3) wystąpi 16 różnych czto- 
nów zawierających Tr GRO) g0)).R(0) [11]. 

Na zakończenie warto podkreślić, że przy meu do przypadku ustalonych baz 
funkcyjnych powyższe wyrażenia na EQ) i E02) znacznie się upraszczają. Zerują się 
bowiem wszystkie pochodne funkcji bazy, co powoduje z kolei znikanie elementów ma- 
cierzy g iS w pierwszym i drugim rzędzie oraz istotną redukcję macierzy fD Di n), 


ZAKONCZENIE 


Rachunek zaburzeń Hartree-Focka jest bardzo dogodną metodą obliczania szeregu 
ważnych magnetycznych i elektrycznych własności atomów i molekuł. W praktycznych 
obliczeniach zmuszeni jesteśmy jednak do zastąpienia dokładnych rozwiązań Hartree- 
Focka przez przybliżone rozwiązania SCF uzyskane w skończenie wymiarowych bazach 
funkcyjnych. Pojawia się wówczas problem zbieżności otrzymywanych wyników do od- 
powiedniej granicy Hartree-Focka. 

W tradycyjnych obliczeniach metodą CHF w ustalonych bazach funkcyjnych zbież- 
ność jest bardzo wolna, co powoduje konieczność stosowania baz funkcyjnych o bardzo 

dużych wymiarach. "pera 
- Wykorzystanie w obliczeniach nie ustalonych, lecz bezpośrednio zależnych od zabu- 
rzenia, baz funkcyjnych istotnie poprawia zbieżność, dając tym samym określone korzyś- 
ci numeryczne. 

Koncepcja zmiennych baz funkcyjnych okazała się bardzo użyteczna w obliczeniach 
magnetycznych własności molekuł [20, 25]. Pierwsze próby wykorzystania tej idei w ob- 
liczeniach własności elektrycznych nie dały w pełni zadowalających rezultatów [21]. 
Przyczyn tego należy szukać w sposobie konstrukcji zmiennych baz zastosowanych 
w tych obliczeniach. 

Ostatnio zaproponowano nową metodę konstrukcji zmiennych baz funkcyjnych, 
odpowiednich do obliczania polaryzowalności dipolowych [22] i kwadr/upolowych [23]. 
Wyniki uzyskiwane w tych obliczeniach pozostają w bardzo dobrej zgodności z danymi 
doświadczalnymi i najdokładniejszymi, znanymi rezultatami metody CHF w bardzo du- 
żych, lecz ustalonych bazach funkcyjnych. Metoda ta będzie omówiona w części II. 


* Przyjęto następującą definicję całki dwuelektronowej 


Kxixk Ix;XD = J f XFL) Xk (2) Ara) XOX) dy, dy, . 
zk 2 
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PE3IOME 


B nauHojń pabote nokasaHo dopMan13M Meroga BO3MylieHuń XapTpu-DoKa MIA 
HOCTOHHHBIX M HEPEMEHHBIX (QyHKUMAHBIX 6a3ncos. PaccMoTpeHo CyYUIECTBEHHYIO 
ripobsiesy CXOĄMMOCTM (KOHBEpPTEHIY1K). Bella noquepKkiiyTa 6onbiiaa cbcPeKTUBHOCTŁ 
6a3M1coB repeMeHHLIX (DYHKIMIi B pacyeTe MATHETMAECKUX M PETREM CBOHCTB 
MHOTO9MEKTPOHHbIX CHCTCM. 


SUMMARY 


In this paper, the Hartree-Fock perturbation theory (CHF) is described both for fixed and variable 
basis sets. For these cases, the problem of ‘convergence is considered. A very high efficiency of the 
variable basis set approach in the calculation of magnetic and electric properties of many-electron 
systems, is stressed. It results from our previous studies. 
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